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Abstract 

Effective astrophysical factors for non-resonant astrophysical nuclear reaction 
are usually calculated with respect to a zero energy limit. In the present work 
that limit is shown to be very disadvantageous compared to the more natural 
effective energy limit. The latter is used in order to modify the thermonuclear 
reaction rate formula in stellar evolution codes so that it takes into account 
both plasma and laboratory screening effects. 
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Nuclear astrophysics experiments are characterized by their continuous effort towards the 
measurement of the zero energy astrophysical factor ( AF) which is used in the thermonuclear 
reaction rate formulas. Recent and past experiments have all attempted to lower 
the energy of the experiment in order to avoid the inaccuracies incurred by the inevitable 
extrapolation to zero energies. That action, however, gives rise to an undesirable electron 
screening effect which enhances low energy data ( |7| and references therein). Moreover, 
experimentalists often fit their AF data with a zero-energy Taylor polynomial which will 
be shown to be less suitable than the more natural effective-energy one. The results of this 
paper should be viewed as an improvement (for some cases) in the theoretical extraction of 
the astrophysical factor and not as a correction of the way experimentalists present their 
results. In the formulas below we follow the formalism and notation of Ref. || . 

The effective astrophysical factor (EAF) S e ff which appears in the non- resonant ther- 
monuclear reaction rate 

7.20 x 1(T 19 N,N 2 2 
ri2= 1 + 5 12 AZ lZ J sTe S ^ (1) 

is defined as 



dE (2) 



and can be calculated by expanding the astrophysical factor S (E) either around zero or 
around the most effective energy of interaction Eq. In the first case the expansion 
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while an expansion around the Gamow peak energy Eq 
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yields( i, 0; 
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where the polynomials are given in Ref. [0. For example, the first two are 

P (n) = 1, P 2 (n) = 12n 2 + 18n + 5 (7) 



In most nuclear astrophysics experiments of non-resonant reaction experimentalist JTT[ 



use a second order truncation of Eq. ([3]) as a fitting formula. They obtain the values 
S (0) , S^ 1 ) (0) , (0) which are then inserted into a second order truncation of Eq. (||) . 
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However, zero energy proximity has serious disadvantages since at such very low energies, 
which are extremely difficult to attain, the experiment is hampered by such effects as beam 
instabilities, impurities, electron screening effects, very small cross sections etc. On the other 
hand the Gamow peak polynomial of Eq. (Q) is hardly ever used although it leads to a more 
accurate calculation of the EAF. In Ref. |H| it was noted that the accuracy obtained by 



the first term of Eq. (H), that is S^jfl — S (E ) , is equivalent to the accuracy achieved with 
knowledge of S (0) and its derivatives. That, of course, is also obvious from Eq. (Q) which 
gives 

OO 1 

S (E ) = £ -SW (0) E n (8) 

n=0 n - 

We can further elaborate on that competition of Eq. versus Eq. by disregarding 
the associated derivatives of the AFs, an approximation which doesn't cause any significant 
error 



12 . In that case we can write 



S^ = S^ + (l + ^)t^(0)ES (9) 

where it is now obvious that S^J incorporates S^j along with an infinite number of correc- 
tive terms. It is again obvious that a simple first order correction to the constant Gamow 
peak AF S (E ) is equivalent to knowing an infinite number of zero energy derivatives. 

On the other hand the crucial region of energies is not close to zero but the Gamow 
window [— A/2 + E Q , E Q + A/2] , with A being as usual the full width at 1/e. Consequently, 
while in the case of S^ff we attempt to lower the energy as close to zero as possible, paying 

(G) 

a great price, when S^j is used the situation improves considerably. First, we only have to 
perform measurements inside the Gamow window, thus avoiding the dangerous zero energy 
region. Then we obtain accurate results by including a limited number of terms while in 
the zero-energy formulas a second order expansion is not adequate as measurements are far 
away from the origin. It is surprising that some experimentalist have been satisfied with 
fitting a second order expansion of the AF given by Eq. ([|) to their data. It is easy to show 
that including some more terms will result in different zero energy AFs. Fortunately the 
results that are used in stellar calculations are often based on independent analyses of the 
experiments (often with theoretical guidance rather than the naive Taylor expansion) 

Another positive aspect of adopting Eq. (||) instead of Eq. @ is the electron screening 
effect, which enhances all very low energy data. The lower the energy, the more pronounced 
the effect which makes Eq. (|5]) seem all the more attractive. In fact, in the above formulas, 
the quantities (0) and (E ) are referring to bare nuclei. However, in the laboratory 
(L) the experimental (ex) values (0) and S^} (E ) are enhanced as they are multiplied 
by the screening enhancement factor (SEF) given by jlT 



h (E) = exp (™|| (10) 
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where n is the Sommerfeld parameter, and is the screening energy obtained via theoretical 

■ e 



models. For nuclei involved in the CNO bi-cycle or more advanced burning stages can 
be given by |14 



U^ = -32.9 Zf /3 Z 2 eV (11) 

where Z\, Z 2 are the atomic numbers of the target and the projectile respectively. Alter- 
natively, if one wants to fully explore the thermal, ionization and relativistic effects on 
screening more elaborate formulas can be used [|T^]. The bare-nucleus AF S (E) is related 
to the corresponding experimental value S ex (E) by the formula: 

S (E) = S ex (E) fl 1 (E) (12) 

Another disadvantage of the zero-energy AF is that we cannot obtain S (0) from the above 
formula, whereas the Gamow-peak AF is readily given by S (E ) = S ex (E ) f^ 1 (Eq). Dis- 
regarding second order derivatives (Eq), which play a minor role and their correction 
(or their inclusion according to Ref. JT^]) would be pointless we have 



Sif f ~ f L (E ) S ex (E ) 



1+ l^ + 2^^ x jEoj Eo + T2^ lnfL{Eo) 



(13) 



Actually, neglecting the screening effect is equivalent to setting (E) = 1. 

The above formula vividly depicts that the value of the EAF used in Eq. ([!]) actually 
depends on the SEF and the model used to describe the screening effect. This of course is 
known empirically to experimentalists who try to correct their data in the laboratory. What 
is derived here for the first time is the analytic dependence of the EAF on the screening mod- 
els selected each time. And there are many, indeed. Another novelty of the above formula is 
that it eliminates all the tedious corrections that are necessary during nuclear astrophysics 
experiments. Moreover, it corrects all nuclear astrophysics experiments performed in the 
past, which have disregarded the screening effect. 

There is another major advantage of Eq. fll"3|) that needs to be underlined. By now it 
is obvious that, if we disregard screening effects, the EAF is approximately equal to either 
S {Eq) or S (0) . Therefore, if the latter is selected by means of Eq. (f|), all temperature effects 
of the EAF are actually neglected, since S (0) will be the same along the whole profile of 
the star. On the contrary, if the former is selected then the EAF takes into account the 
region of the star where the reaction takes place via the temperature dependence of the 
most effective energy of interaction Eq (T) . For example, if we adopt Eq. (|j) in our solar 
codes using the bare-nucleus value || S (0) = 5.32 for the He 3 (He 3 , 2p) He 4 reaction then 
that will practically be the value of the EAF at the center of the sun, throughout its energy 
production core as well as in any other stellar interior no matter how hot it is. That is indeed 
an unnatural result. The only way Eq. (Q) could counterbalance the temperature dependence 
of the leading term S (Eq) of Eq. (Q) is to include an infinite number of temperature 
dependent terms, which is of course futile since that can be accomplished by merely setting 
S e ff — S (Eq) in the reaction rate. 

To gain an idea of the screening corrections to the EAF let us consider the break-up 
reaction He 3 (He 3 , 2p) He 4 whose SEF according to a recent model [II is fi (E ) = 1.23. 
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The screening correction inside the brackets in Eq. ([Op is roughly (15/12r) In (E ) = 
4 x 10 -3 while the first significant term is (5/12r) = 7.7 x 10~ 3 . Obviously, the screening 
correction is comparable to the first corrective term which is usually retained in the formula 
of the EAF. Of course they are both negligible and the only significant correction in Eq. 
(|j~3|) is that of the SEF multiplied by S ex (E ) . As for reactions involved in more advanced 
burning stages than the pp one we can consider the most important reaction of the CNO 
bi-cycle namely N 1A (p, 7) O 15 . For that reaction the two previously mentioned corrections 
are roughly the same (~ 6. x 10~ 3 ). Note that including any of the corrective terms inside 
the brackets would be meaningless for another reason. As we observed their contribution is 
of order 1% while the experimental error of the leading terms 5* (Eq) (or even S (0) when 
Eq. (f|) is adopted) is much larger 0. 

We can now safely argue that all screening effects on the EAF can be taken into account 
by multiplying the (uncorrected) EAF with the SEF given by fi (Eq) . Therefore the re- 
action rate itself is now multiplied by a SEF which is the combination of laboratory and 
plasma screening effects. The reaction rate is now written 



7.20 x IP' 19 N,N 2 2 (G) 
n2= 1 + 6 12 AZ lZ J SLTe W (14) 



where S^jJ is merely a single experimental measurement, that is S^J ~ S ex (E ) , and the 
combined SEF f$L is given by 



fsL (E ) = exp 



Tin (E ) 

En 



{u! - ul 



(15) 



Note that the values of the screening energies (u^, U^j in the above formula are absolute 
ones. At first sight it seems that Eq. (|I4"D and Eq. flTSP indicate that r 12 7^ r 21 , since 
(Z\, Z 2 ) 7^ Z7g (Z 2 , Zi) , an error also committed in the plasma SEF derived in Ref. ||16 



This is not the case here. The combined SEF given by fsL must always be coupled with 
the appropriate experimental value S ex (Eq) so that it is always r\ 2 = r 2 \ . 



According to Eq. (|l5j ), plasma and laboratory screening have opposite effects on the 
reaction rate. The theoretical value of has already been discussed. On the other hand 
the value of the plasma screening energy can be derived using various theoretical models 
such as Salpeter's [|I7| according to which = Z\Z 2 e 2 r^ where the Debye radius, or 
even Shaviv's |1§ according to which Jjf =l.hZiZ 2 e 2 r]^ . 

Regarding these two plasma screening prescriptions a comment is imperative. In a recent 
work Shaviv's prescriptions has been argued against. However, the relation between 



Salpeter's approach and Shaviv's is as it should be. Actually they provide upper and lower 
limits for the screening energies, as has often been done in laboratory screening models fl4" 



where a sudden (Salpeter's) and an adiabatic (Shaviv's) limit is defined. It is important to 
realize that the most accurate way of deriving a SEF is to provide a sudden and an adiabatic 
limit with a minimum discrepancy between them. The two formulas in question fulfill that 
requirement and constitute the most reliable screening constraints at the moment. 

In conclusion, the competition between Eq. (|| and Eq. (f|) is an unequal one. Just a 
single term of the former yields a more accurate EAF than a great number of terms of the 
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latter. Moreover, Eq. (|6|) takes into account temperature effects in a consistent and accurate 
way, whereas Eq. (Q) practically disregards them. Finally, using Eq. we were able to 
modify the reaction rate formula so that it automatically takes into account the laboratory 
electron screening enhancement. The combination of the plasma and laboratory screening 
led to the derivation of a combined screening factor which allows theoreticians to have full 
control over the theoretical models used for the description of the electron screening effect. 
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